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Abstract. We present a new algebraic extension of the classical MacMahon Master Theo- 
rem. The basis of our extension is the Koszul duality for non-quadratic algebras defined by 
Berger. Combinatorial implications are also discussed. 



Introduction 

The MacMahon Master Theorem is a fundamental result in Enumerative Combinatorics, 
with a number of applications to partition theory, binomial identities and q-series. Ever since 
Percy MacMahon introduced it in |MM| . the theorem was understood largely in a combinato- 
rial and analytic context (see |CF1 IGJj l. In an important development, a quantum analogue 
was discovered by Garoufalidis, Le and Zeilberger jGLZj . sharply extending a basic partially 
commutative version of Cartier and Foata |CF| . Since then, several other algebraic extensions 
have been introduced |FH21 IRT] . but their nature remained mysterious. 

Most recently, Hai and Lorenz showed that the quantum analogue can be obtained from 
the Koszul duality |HL| . While |HL| does not rederive the results of other extensions (which 
appeared independently), it goes a long way to explain them, in particular the fact that all 
these algebraic extensions are based on quadratic algebras with Koszul property. It is a natural 
question whether further (non-quadratic) extensions can be obtained. 

In this paper we present an extension of the MacMahon Master Theorem of a new type. The 
construction is based on special non-commutative algebras B m ^ studied by Berger in [B]. These 
algebras have m generators and homogenous relations of degree k coming from the exterior 
algebra. Algebras B m ,k are fundamental examples of the generalized Koszul property for non- 
quadratic algebras defined in [B]. Our proof is uses the Hai-Lorenz approach and Berger's 
explicit calculations of the Koszul resolutions in this case. The result is a curious algebraic 
identity which can be stated in a combinatorial language. 

In this note, we first recall the MacMahon Master Theorem and state the main result in 
Section^ We then give a combinatorial interpretation in Sectional Two natural examples and 
enumerative applications are considered in sections |3| El while the proofs are postponed until 
Section [5] We conclude with final remarks and open problems (Section [HJ). 

In conclusion, let us say a few words about the fundamental difference between an algebraic 
and a combinatorial approach to the Master Theorem. A combinatorial approach starts by 
translating the algebraic identity as an equation for two sums of words in a formal language. 
This equation is then proved by an explicit bijection. The algebraic approach starts by showing 
that both sides of the algebraic identity are invariant under a large group of symmetries. This 
action is then utilized to obtain a concise and insightful proof of the result. The main difference 
is that in the combinatorial approach the bijection itself is of value, and its properties often 
lead to new results. We refer to jCFI IFZI IKPj for examples of this phenomenon in the context 
of the MacMahon Master Theorem. It would be interesting to obtain a combinatorial proof of 
our main result (Theorem |2J) and further extend it in these directions. 
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1. Main results 

We begin by stating the MacMahon Master Theorem in the classical form: 

Theorem 1. (MacMahon Master Theorem) Let A = (a.;j) mxm be a complex matrix, and let 
Xx, . . . , x m be a set of variables. Denote by G(k±, . . . , k m ) the coefficient of x^ 1 ■ ■ ■ x^ in 



(1) Y[{anxi + . . . + a im x 



ki 



t=l 



en 



(4) ^(-lr^^.-.^^o, 



Let t±, . . . , t m be another set of variables, and T = (Sijti) mxrn . Th 
(2) (fei E m) ^--->U^---4- = ^T^, 

where the summation is over all nonnegative integer vectors (ki, ■ ■ ■ , k m ). 
By taking t\ = . . . = t m = 1 we get 

(k u ...,k m ) y ' 

whenever both sides of the equation are well defined. It is easy to see that (01 and J2J hold 
when all Oij are formal commutative variables. Moreover, replacing a.^ in © with a^ti shows 
that is equivalent to @ , so in the literature the MacMahon Master Theorem is often stated 
in this form. 

Let us move now into a non-commutative setting. Fix an integer k, 2 < k < m. Let B = B m .k 
be an algebra with basis X — {x\, . . . , x m } and the following defining relations: 

gnO) 

for all 1 < ii < . . . < i}~ < m. When k = 2 we get the usual symmetric algebra C[xi, . . . , x m \ 
Define the set A(m, k) of admissible sequences i = {i\ . . . if ), where i\, . . . ,ij> € {1, . . . , m}, and 
such that no k subsequent indices are strictly decreasing. We call monomials 

Mi(X) = Xi x ■ ■ ■ Xi e , where i = . . . ,ii) £ A(m, k) , 

the admissible monomials. It is known (see p. 723]) and easy to see that the set of admissible 
monomials Mi{X) is a linear basis in B (use a Grobner basis argument with the lexicographic 
ordering). When k = 2 we get the usual basis in C[xi, . . . , x m ] with indices . . . , if) € A(m, 2) 
weakly increasing. 

Now fix a matrix A — (dij), £ C, and let 

Mi(AX) = y il ■ ■ ■ y l( , where yi = anxi + . . . + a. lm x m , for all 1 < i < m. 

The product Mi(AX) is an element in algebra B and is a generalization of the product For 
an admissible sequence i = . . . , it) £ A(m, k), denote by G(i) = G{i\, . . . ,ii) the coefficient 
of Mi(X) in Mi(AX). We set G(i) — 1 for an empty sequence i = 0. 

Theorem 2. Let A = (aij) mxra be a complex matrix, let t\, . . . ,t m be another set of variables, 
which commute with Xi and each other, and let T — (5ijti) mxm . Then 

(5) G{i u ...,u)t il ---t il ■ ]T {-IT^cATA) = 1, 

(ii,...,i«)eA(m,fc) r = 0,l modfc 

where a(r) = r — (rmodk), and polynomials Ci(TA) £ C[ti, . . . ,t m ] are defined as coefficients 
of the characteristic polynomial: 

m 

(6) det(A7 -TA) = J2 c r (TA) \ m - r . 

r=0 
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Note that when k = 2 the power a(r) is always even, and r takes all nonnegative values < to. 
Therefore, the second sum in © is equal to det(7 — TA), as in |3J). Similarly, each term in the 
first summation in (JjjJ is equal to the corresponding term in J2J. This implies that Theorem [21 
is an extension of the MacMahon Master Theorem ^ 



2. COMBINATORIAL REFORMULATION 

In this section we restate our main result, Theorem |3 in the language of words. We assume 
that ay , 1 < i,j ' < m, are commutative variables. We begin with some combinatorial definitions. 

Denote [to] = {1, . . . , to} and let i = . . . , i{) G [m] . Recall that A(m, fc) denotes the set 
of admissible sequences i. Denote by A(m, fc, £) the set of sequences of length £. Let 

inv(i) = \{(s,t) : i s > i t , s < t}\ 

be the (usual) number of inversion in sequence i. Define the smallest decreasing k-sequence in i 
to be a subsequence i s > i s -\-i > . . . > i s +fc_i, where s is the smallest possible. 

We say that two sequences i, i' G [m] e differ by a k-reversion, write i — ► i', if i differs from i' by 
a permutation of elements of the smallest decreasing fc-sequence in i'. Of course, inv(i') > inv(i) 
in this case. For example, (4,3,2,6,1) — > (4,6,3,2,1), where k — 3, to — 6, £ — 5, and 
(6, 3, 2) is the smallest decreasing 3-sequence in (4, 6, 3, 2, 1) 6 [6] 5 . Here inv(4, 3, 2, 6, 1) = 7 < 
inv(4,6,3,2,l) = 9. 

By a reversion path from i to j, where i, j G [to]^, we call a sequence of fc-reversions: 

7 : i - i' - i" - ... - j 

We write i — >j j in this case, and let I7I denotes length (number of fc-reversions) in path 7. Since 
at every fc-reversion the number of inversion increases, we clearly have inv(j) — inv(i) > |7|. 
Define integer coefficients 

(7) CiJ = (-l) inv (j)- inv ( i HTl , 

where the summation is over all reversion paths 7 from i to j. Note that when k = 2, we have 
Cj.j G {0, 1} depending on whether j >- i in the weak Bruhat order. Let us consider the products 

a ij = a iiji ' ' ' a iiji i 

where i = (%%, i e ), j = (ji, . . .,j t ) G [m] e . 

A partial permutation lo of [to] is defined as a pair J C [m] and the usual permutation 
of J (which we also denote by to). By inv(w) denote the (usual) number of inversion of the 
permutations of elements in J. Denote by Y, m (r) the set of partial permutations of [m], such 
that I J = r. For a partial permutation lu G S m (r) define 

where J = {ji, . . . ,j r }- Set a w = 1 for an empty partial permutation (when J = 0). 
Corollary 3. For k < to, in notation above we have: 



(8) 



1=1 i£A(m,M) jeH' 



E (-i) Q(r) E (-!) invM «- =1, 

r = 0, 1 modfc u)£S m (r) 

where a(r) = r — (r mod k). 

Let us emphasize that the first term in @ is an infinite sum. In contrast with the case 
fc = 2, it no longer has only positive coefficients. On the other hand, the second term is a finite 
summation, and the sign pattern is completely prescribed. 

In fact, the corollary is equivalent to Theorem[21 This follows from the proof of the corollary 
given in Section [5j 
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3. First example 

Let A = (dij) = I be the identity matrix. Then Mi(AX) = Mi(X) and G(i) = 1, for all 
i G [m] . Therefore, the first term in j£| is equal to the weighted sum over all admissible 
sequences, which we denote by F mj k ■ 

Fjn^k (tl , ■ ■ ■ , ^m) — ^ ^ tii ' ' ' tig • 

(ii,...,ii)eA(rn,k,£) 

For the second term, note that 

det(AJ - TA) = (A -<!)••• (A - t m ). 

Therefore, the coefficients of the characteristic polynomial defined in JSJ are equal to the el- 
ementary symmetric polynomials: Cj = (— l)'ei(ti, . . . , ti). Now equation (JSJ gives a closed 
formula for the infinite summation F m ^ ■ 

(9) F m , fc = ' 



1 — ei + efc - efc+i + e 2 k - e 2 k+i + ■■ ■ 

where indices in the denominator are < m. In particular, this implies that F m ^ is a symmetric 
function (see SectionElfor further remarks). Of course, when k = 2, equation (j^J is the classical 
symmetric functions identity connecting complete and elementary symmetric functions: 

do) i +hl+h2+h3+ ... = ^^j.,,^)^ - 

We refer to |Mal IS] for the background, proofs and other results on symmetric functions. In 
the opposite extreme, when k = m, we have: 

1 



Fm,m(tl,, " ,g " l-(*l + ... + U + *l---*m' 

This equation can be shown directly by a simple inductive argument. 

Consider now the enumerative applications of our equations. Denote by L m ^(l) = |A(m, k, £)\ 
the number of admissible sequences. Taking t\ = . . . = t m = t in © we obtain a generating 
function for the number of admissible sequences: 

(12) l + JW^^l"^^-^^ 1 

While the exact shape of (|12fl may seem surprising, the fact that it is rational function is 
not. To see this, consider a graph T m ,k on A(m, k, m) with (directed) edges of the type: 

(ii,i 2 ,...,i m ) — ► (i2,h, ■ ■ ■ ,i m , j) ■ 

Clearly, L m ^{m) = m l — (™), since (i\> i 2 > ■ ■ ■ > i m ) are the only non-admissible sequences 
in this case. Observe that admissible sequences of length I correspond to paths of length 
(£ — m + 1) in T m . Using the transfer matrix method (see jSJ Chapter 4]) we obtain a rational 
generating function for L m ^{m), for any fixed m > k. Note also that the transfer matrix 
method gives denominator of degree |r m | ~ m e , much larger than the degree in (|12|l . where it 
is < TO. 



4. Second example 

Let k — m and suppose A — (ay), where <2y = 1 for all i,j G [to]. In this special case we 
have: 

Mi(AX) = ( Xl + ... + x m ) 1 , 
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for all i G [to] . Setting t% = . .. = t m = t in ((SJ, we obtain the first term: 

oo 

Y, G(i)t l = 1 + Y N ™(£)*> 

i£A(m,M) 1=1 

where N m (£) is the sum of the coefficients of all admissible monomials in [x± + . . . + x m ) 1 . 

Now, the module condition for the second term in (0 are very restrictive: only r = 0, r = 1, 
and r = m are allowed. There is a unique empty partial permutation corresponding to r = 
case, and exactly to partial permutation corresponding to r = 1 case (all consisting of single 
elements). Theorem[2]now gives: 

00 1 I 

jr[ m{ ) |E m (0)| - |E m (l)|t + det(A)t™ 1-mf 

We conclude that N m {£) = m e , which may seem quite surprising given that the number L m _ m (£) 
of admissible sequences in A(m, to, £) is much smaller than m l . To see this, recall the generating 
function from the first example: 

oo ^ 

da) i + ew^= 1 _ mt + tm , 

1=1 

which is a special case of 11110 . From ()13f) it follows easily that L mm (£) — + o(l)) as 
£ —¥ oo, for some j3 < to. Of course, this can also be seen by a direct argument. 

Now, to explain the identity N m {£) = m e , consider what happens when the relations l@J in 
algebra B m ,m are applied to the to terms in (x\ + . . . + x rn ) £ . We have: 

(14) ^=^(-1)^^, 

j' 

where the sum is over m! — 1 permutations j' of the smallest decreasing subsequence in j and 
the sign is the sign of these permutations. Clearly, the total sum of coefficients of terms on the 
right-hand side of (|14fl is equal to 1. Therefore, the sum of coefficients of terms is unchanged 
under the algebra relations, and thus equal to m , the initial sum of all coefficients. 

5. Proof of results 

Proof of Corollary\Q First, note that we used equation (JSJ with t\ = . . . = t rn = 1. The second 
term in (jSJ can be easily seen to be equal to the second product in J5J|. This follows from the 
usual expansion of determinant: 



det(XI-A)=J2 Y A" 1 r (-I) 1 " - 



\mv(uj) 

r=0 uj£S m (r) 



as in JBJ . Now, to see that the first term in JSJ is equal to that in (J5J we first observe that 
(15) Xj 1 ■ ■ ■ Xj t = Cij x% t ■ ■ ■ Xj e , 

i=(ii,...,if)eA(m,fc,£) 

for every j = (ji, . . . ,je) G [m\ l . This follows by induction. Start with a product Xj t ■ ■ -Xj e 
corresponding to j and apply the relations in algebra B m ,fc to the smallest decreasing fc-sequence. 
We obtain products Xj^ ■ ■ ■ Xj^ corresponding to k- reversions xy — ► Xj . Now apply this relation 
again, and repeat this until only admissible sequences are obtained. Then each reversion path 
i — j corresponds to a term in the summation. After checking that the signs in @ correspond 
to those given by J3J, we obtain equation l|15|) . 

Note now that G(ii, . . . defined as the coefficient of Mi 1: , ^i t {X) in M(AX), is equal to 
the summation of products ay with coefficients cy. Checking the terms of the summation, we 
obtain the result. □ 
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Proof of Theorem^ Before we move to the general case, let us first recall the algebraic proof 
of the MacMahon Master Theorem (see |GLZI iHLj ). We rewrite equation © as follows: 

oo - 

Y Tr (S i A) = - l T - , 
det(Z-A)' 

where V denote the vector space of variables X\, . . . , x m . Here the left-hand side is equal to the 
left-hand side of © by definition of the symmetric power S (A). Now, there is a natural action 
of GL m (C) on both sides. Since matrices with distinct eigenvalues are dense in all matrices, 
it suffice to show the result for diagonal matrices. In this case equation J3J is equivalent to 
equation J2J) for A = I. Finally, the latter can be written as (|10f> which can be proved by a 
straightforward calculation (see e.g. Ma, S ). 

Let us note here that S e V is a homogeneous component of the natural action of GL(V) 
on the symmetric algebra A = C[x%, . . . ,x m ]. Then the left-hand side of is exactly the 
character of the GL(V)-&ction in A. Now Theorem 1 follows immediately from here. 

In the general case k > 2 we can proceed in a similar way and reduce the theorem to the case 
of diagonal matrices, which in turn reduces to the case A — I. As shown in Section [21 this can 
be written as © ■ Unfortunately, the above mentioned straightforward calculation is no longer 
possible in this case (cf. Section IB3|) . Following the idea in |HL| . we take a different approach. 

We use the notion of generalized Koszulity defined in [E]. By Theorem 3.13 of the 
algebras B defined in Section ^ are generalized Koszul. We can now proceed with the proof of 
the theorem. 

First, note that it is sufficient to prove the theorem in the case when A is a nondegenerate 
matrix, since such matrices are dense in the space of all matrices. Thus, we may assume that 
A G GL m (C) = GL(V), where V (as before) is the space of generators x\, . . . , x m . 

Since the relations of the algebra B are GL(F)-invariant, the homogeneous components of 
this algebra are representations of GL(V). Therefore, the characters of these representations 
may be computed from the Koszul resolution (see |F3 IPP| for the background). It follows 
from |2] that 

To4 j+£ (C,C) = A kj+E V, 

sitting in degree kj + e, where £ = 0,1, and where A V denote the exterior powers. Therefore, 
the computation of the character via the Koszul resolution yields 

^2Tr(A\ B[p] )t? = (j2 E (-l) £ Tr(A| AfcJ+ey )^ 

p>0 \j>0 £=0,1 

where B[p] denotes the p- homogenous component. Recall that the first factor in (JSJ is exactly 
the character of the GL(F)-action in B. The above equation then implies the theorem. □ 




6. Final remarks 

6.1. There are other examples of algebras with k- homogeneous relations which satisfy general- 
ized Koszul property 0. For example, one can omit all signs in Q and obtain a result similar 
to Theorem |2 where the second summation in |j^J has to be modified by replacing exterior 
powers with symmetric powers. By analogy with l|12f) we get the following generating function 
formula: 

( il ,...,^ (ro ,M) l-hx + h k - h k+1 + h 2k h 2k+1 +... 

where A'(m, k, I) denote the set of sequences (ii, . . . , it ) 6 [m] e without k subsequent elements 
which are weakly decreasing. When k = 2 we yet again get up to a substitution i, <— (—if). 
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For k > 2, both the left-hand side and the denominator in the right-hand side are infinite series, 
so the result is somewhat less natural from combinatorial point of view. 

6.2. Following |OLZ| (see also |FHlllHLllKP| 1. there is a natural quantum analogue of Theo- 
rem The algebra relations J3J have to be replaced with 

where q € C, q ^ 0. This algebra is a flat deformation of B m ^, and its generalized Koszulity 
follows from the Grobner basis argument in Berger's proof for B m ^ 0- Another difference 
is in the equation © which will now contain g-minors. This generalization of Theorem |2) is 
equivariant under U q {GL{V)) and can be proved along the same lines. 

One can similarly find fc-homogenous extensions of the Cartier-Foata Theorem and other 
recently considered generalizations of the MacMahon Master Theorem (see |FH2I IKP") 1. 

6.3. Although equations JHJ and (|16|l do not seem to be available in the literature, there are at 
least two ways in which they can be derived from the much more general known combinatorial 
results. On one hand, the equations follow from the Maximal String Decomposition Theorem 
(see |GJI §4.2]), and on the other hand from the Cluster Theorem (see |GJI §2.8] and |5) 
Exc. 4.14]). It would be interesting to find a direct bijective proof of the equations. 

As we mentioned in Section 13) equation @ and the eigenvalue argument gives an alternative 
proof of Theorem [21 a proof which does not use Koszulity. Moreover, one can reverse the 
argument in the proof and derive Berger's result on the Hilbert series of B m j. from the theorem. 
In other words, one can now prove Berger's result combinatorially, from the Cluster or MSD 
theorems. 

Finally, let us mention that both Cluster and MSD theorems work in the non-commutative 
setting. In particular, the second part of the MSD Theorem (as stated in |G.I| 1 gives a natural 
g-analogue according to the number of inversions in a sequence. We are curious to see the 
relationship between this q-analogue and the quantum analogue described above. 



6.4. To see that F m ^ is symmetric simply observe that 

E 



Fm,k (*£l ; • • • 7 ) — / Su\<El : • • • ; <£m) ; 



where the summation is over all finite diagrams v which start at the upper left corner and move 
down and to the right, and such that all columns have size smaller than k. Such v are not 
to be confused with the rim hooks, which are reflections of these diagrams (see e.g. |Mal 15)1. 
Although such v are not skew Young diagrams, Schur polynomials s„ can be defined for any 
finite arrangement of squares in N 2 . In this case s v is a summation over all sequences with the 
descent set corresponding to v. Their symmetry follows from the usual considerations. 

6.5. Even though the equation Ijl2(l is a simple evaluation of JjJJ, the closest we could find in 
the literature is the following formula: 

^ t n ( x k x k+1 

1 + > en -r = 1 - X 



^ n n\ V ' kl (fc + 1)! 

n— 1 x x 7 

where e„ is the number of permutations in S n with no decreasing runs of length k (see e.g. |DB| 
and jC.II Exc. 4.2.8]). The connection between enumeration of permutations and sequences 
with forbidden subwords is also known (see |OTI §4.2] for details and further references). 

6.6. As we mentioned in the introduction, the most famous application of the MacMahon 
Master Theorem is to binomial identities |GT| . It would be interesting to find applications of 
Theorem|2in this direction. 
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